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(N Abstract 

C^) Copper and aluminium sub-micron films are structured as columnar arrays of grains traversing the entire film, and hence are virtual 

O^l two-dimensional polycrystalline solids. A closed-form solution for the two-dimensional plastic flow is obtained assuming that 

S— i grain boundary sliding is the dominating mechanism for deformation. The theoretical results obtained here are in close agreement 



cd 



with experimental data appeared in the literature, which suggests further work on the subject. 
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Tensile and compressive stresses of thermal origin can reach 
hundreds of MPa in the small-scale structures of a complex 
microelectronic device, either on service or in the fabrication 
process. Besides, the strong functional interdependence of the 
numerous constituents of the system leaves no tolerance for 
local strain failures. The structural dimensions of current mi- 
crosystems are in the length scale of 100 nm or less. Hence, 
reliable design ensuring no risk of stress-induced catastrophic 
defect demands precise knowledge of the mechanical properties 
of sub-micron metallic films, together with a good understand- 
ing of the response of the film to strongly varying stresses. This 
has opened a new field of research in materials science because 
the mechanical behaviour of thin films has proven to be very 
different from that of their bulk counterparts. 

Normal-view and cross-sectional micrography has shown 
that polycrystalline sub-micron films have grain sizes larger or 
of the order of the film thickness, with grains traversing the en- 
tire film |Q~1[2)- Typical films are structured as columnar arrays 
of grains which in the plane of the film exhibit random equiaxed 
shapes and crystallographic orientations. However, partial crys- 
tallographic ordering is detected in the direction normal to the 
film plane. The grains of copper sub-micron films, either self- 
standing or on a supporting substrate, show strong tendency to 
orient the (111) crystal direction along the normal to the film 
surfaces (UG). 

A sub-micron metallic film of this kind configures a physi- 
cal realization of a two-dimensional polycrystalline solid. The 
interest of this goes beyond the technical applications because 
provides means for the laboratory testing of theoretical mod- 
els for plastic flow in a reduced dimensionality. Techniques for 
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the fabrication and mechanical testing of self standing beams 
of sub-micron thickness have been developed in recent years 
0. Besides this, a neat solution for the problem of drawing 
out the intrinsic film mechanical properties when testing thin 
films on substrates has appeared quite recently in the literature. 
In it the true stresses exerted in both principal directions of a 
metallic film, deposited on a compliant polymer substrate, are 
accurately measured by diffracting in situ the X-rays coming 
from a synchrotron source 0. The diffraction pattern reveals 
the distortions of the crystal cells of the metallic film, which 
are proportional to the true stresses operating on the film. The 
stabilising effect of the substrate permits relatively large plas- 
tic tensile and compressive deformation of films as thin as 80 
nm without the risk of buckling. Because of the mismatch of the 
Poisson ratios of the metallic film and the substrate, any strain s 
in the longitudinal direction is associated not only to a stress cr L 
in that direction, but also to a stress o~t in the transversal one. 
This is not a serious shortcoming because the method gives both 
stresses, <tl and <Tt, and what matters is the difference cr L - crj 
between them, as will become clear in what follows. 

Experiments on copper [2| and aluminium [3j micron and 
sub-micron films show that the plastic deformation starts with 
a strong but brief hardening period, going to not beyond 0.5% 
strain, followed by a plateau or a slowly softening regime ex- 
tending up to failure. Simulations by von Blanckenhagen et 
al. [4] of the response to an applied stress of the dislocation 
sources in a representative columnar grain, and the subsequent 
dislocation dynamics, explains very well the initial strengthen- 
ing plastic deformation, fitting both the measured stress-strain 
curves and dislocation densities. However, for strains exeeding 
0.5% the simulation gives too large stresses and does not ex- 
plain the plateau, suggesting a crossover with another deforma- 
tion mechanism. We provide here supporting evidence that the 
latter higher strain mechanism is two-dimensional grain slid- 
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ing. 

If we assume that grain boundary sliding is the dominant 
mechanism for the plastic deformation of the film, we can 
reduce the general theoretical formalism of Ref. to two- 
dimensions for analysing the data. The basic hypothesis of the 
model is that adjacent grains can slide past each other over long 
distances by effect of the shear stresses actuating in their shared 
boundaries, accommodating at the same time their shapes by 
internal mechanisms to prevent voids at the interfaces and pre- 
serve matter continuity. The sliding-induced stress fields as- 
sociated to grain shape accommodation are of little relevance 
because are assumed much weaker than the shear stress caus- 
ing grain sliding. The latter is linearly related with the relative 
speed |Av| between the adjacent grains and has a threshold t c , 
below which no elementary sliding process can occur |6|. 

Fig. [T]represents a pair of adjacent two-dimensional grains. 
The x' and y' axes of the local frame of reference (x'y'), asso- 
ciated to the unitary vectors f and j', are normal and parallel to 
the grain boundary shared by the grains, respectively. The main 
frame of reference (xy), with unitary vectors i and j, has its axes 
in the principal directions of the stress tensor. Denoting cr,-y, 
/',/ = x',y', the components of the stress tensor in the (x'y') 
frame of reference, the relative velocity Av of the two grains is 



Av = 



(<2(<x v y -s't c )]' if cry,,, > t c 



0, 



otherwise, 



(1) 



where Q is a proportionality coefficient and s' = cr v y/|oy v ,| 
is the sign of the shear stress cry,,,. The term s't c ensures that 
Av = when \<T x > y >\ - t c . This expresion for Av has proven to 
hold with great accuracy for several aluminium, titanium and 
magnesium alloys J7]. The coefficient Q must not depend on 
either the shear stresses or the orientation of the grain boundary, 
therefore its dependence on the normal stresses is only via the 
hydrostatic pressure invariant p - -(cr x , x , + cr v y)/2 |8, 9|. 
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Figure I: Schematic representation of two adjacent grains showing the local 
(x'y') and main (xy) frames of reference, the corresponding unitary vectors and 
the relative velocity Av. 



Replacing <x A .y = -(°~x _ 0y) sin cos 9, where cr x and <r y 
are the principal stresses and 6 is the rotation angle defined in 
Fig.fT] and j' = - sinfli + cosf?j, it gives 



Av* = Q[-(cr x - a- y )sm8cos0 + s'T c ](-sin0i + cosf?j), (2) 



which presumes that 



sin(26>) > 



2t c 



= sin(2f? e ), 



(3) 



and otherwise Av = 0. A grain boundary whose normal i' sub- 
tends with the principal direction x or y an angle smaller than 
9 C is not able to slide because the in-plane shear stress is below 
T c , no matter how strong the external forces may be. 

To link these equations with the velocity field v(x,y) of 
the plastically flowing film, consider two points at (x,y) and 
(x + 6x,y). The macroscopically small segment 6x intersects a 
large number n of grain boundaries, and then dx - nd, where d 
is the mean grain size. The relative velocity between the start- 
ing and final points of the segment 6x is the sum of the n rela- 
tive velocities between the consecutive grains it passes through. 
Thus 



v*(x + cue, y) — v(x, y) 
Sx 



A-=l 



(4) 



where k numbers the successive grain boundaries intersecting 
Sx. But for the factor d in the denominator, the right-hand-side 
of this equation defines an average. Consequently, in the proper 
limit, 

dv 1 

— = -7<Av) ; , i = x,y or x t = x,y, (5) 

ox t d 

where symbol (...), means the average over all boundary ori- 
entations Q compatible with x, > 0. With the latter restriction 
and recalling condition pi) this means 



iff 
dxi 



nd J D . 



d6 Av, i — x,y or x; — x,y , (6) 



where the integration domains D x and D y are those shown in 
Fig. [2] Solving one obtains that 



dv x 


dv y 2Qt c 


t cos 3 6 C - sin 3 6 C 


dx 


dy nd 


\ 3 sin 6 C cos 8 C 


dv x 


dv y 
dx 




dy 





cos 6 r + sin 9, 



where 



= - arcsin 



2t> 



■cr,\ 



(7) 



(8) 



These equations assume that the mean grain size d is con- 
served, and yield a complete set of equations for the velocity 
field v of the material medium in terms of the stresses. From 
them 



Vx 



A 
A 



(9) 




Figure 2: The shaded areas represent the integration domains D x (a) and D v 
(b) for the variable in integrals JSl. 



which means that the plastic flow is laminar and conserves the 
film specific area A. The strain rate tensor sy = (dvj/dxj + 
dvj/dxj)/2 follows directly from Eqs. (FTl), explicitly 

2Qt c /cos 3 6 c - sin 3 6 C \ 

Sxx = -£yy = —7- i ■ 7, COS 6 C + Sin 6U , 

nd \ 3sin0 c cos0 c / (10) 

s xy = 0, (<r x > cr y ) 

for the biaxial deformation with principal axes x and y. 

Eqs. ( 10 1 and ((HI) show two very important properties, spe- 



cific to the two-dimensional plastic flow by grain sliding. They 
are: 

(i) The stresses depend only on the strain rates, and not on the 
strains. Hence no strain hardening or softening of the material 
is expected in a bidimensional grain sliding plastic regime, 
(ii) The ratio between the true strain rates and the coefficient Q 
depends on the principal stresses only through their difference 
0~ x - cr y . As explained before, the coefficient Q is expected to 
depend on the pressure invariant p = -(cr x + <x v )/2, however 
this dependence seems to be rather weak at room temperature, 
as revealed by the experimental data shown next. 

Fig.[3]shows the data of Gruber et al. [2] on a copper film of 
67 nm thickness and mean grain size d - 104 nm (measured in 
the film plane), deposited on a polyimide substrate and strained 
at constant strain rate. The film incorporates a capping tantalum 
layer of 9 nm thick. The principal stresses 07, and cr T , which 
can be assimilated to cr x and cr v in Eqs. ( 10 1 and d8jl, were mea- 
sured by the synchrotron X-ray diffration method. The graph 
shows that both cr L and cr T vary with the strain e, which is 
assimilable to e xx in the previous equations. The longitudi- 
nal stress rapidly increases to a maximum <x L = 846 MPa at 
e - 0.027 and then decreases monotonically to ctl = 515 MPa 
when the tensile test ends at s - 0.066. The transversal stress 
crj- displays a similar behaviour, with maximum cry - 404 MPa 
at s = 0.027 and cr T = 140 MPa at s = 0.066. Both stresses 
vary a few hundreds MPa and are far from constant. 

Instead, the difference <ji - crj (white circles) increases 
rapidly in the interval < e < 0.015 and then remains nearly 
constant up to the end of the tensile test at e = 0.066. Recall- 
ing Eqs. (10 1 and (fHJ, this evidences that the tensile test was 
carried out at a constant strain rate. To be precise, the cross- 
head displacement of the mechanical testing machine was in- 
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Figure 3: Data of Gruber et al. |2| for a copper film of 67 nm thick on a 
polyimide substrate, strained at constant strain rate. The film has a passivating 
tantalum layer of 9 nm thick. Crosses represent the original data and circles the 
difference ctl-ctj between the two principal stresses. Within the experimental 
uncertainty, o - /, - crj in the tensile region is constant for £ > 0.01. 



creased in small steps because the registration of a diffraction 
picture takes one or two minutes. Also, the flat region of the 
curve <j L - cr T vs. s has a slight negative slope, which can be 
explained because what was kept constant in the experiment 



was the engineering strain rate, while s xx in Eq. ( 10 1 is the true 
strain rate. When the engineering strain rate is maintained con- 
stant, the true strain rate diminishes with strain, and the same 
does cr x - <x v , as given by Eq. ( 10 1. It is usually argued that 



plastic deformation is driven by the deviatoric principal stresses 
<x, + p, where the hydrostatic pressure p is the negative average 
of the normal stresses. In two dimensions the deviatoric stress 
is (cr v - cr v )/2. In the insets of the figures we disregarded the 
factor 1/2 and called the deviatoric stress to <x L - crj- 

Upon unloading, Fig. [3] shows that <x L - cr T varies linearly 
with e, with a slope close to the one displayed by the elastic part 
of the loading curve. On compression, cr^-crj converges rather 
slowly to a compressive net stress 419 MPa, which is a quantity 
amazingly close to the tensile strain <x t - crj at the same strain 
in the tensile plateau of the experimental curve. This shows that 
the material undergoes no final strain strengthening or weaken- 
ing, that is, no Bauschinger effect seems to occur. 

Figs.[4]and[5]displays the data given by similar experiments 
with copper films of different thicknesses and mean grain sizes, 
and with one or two passivating tantalum layers of 8 nm thick. 
As in Fig. [5] and by the same reason, in both cases the plateau 
exhibits a small negative slope. Upon compression, \<j L - cr T \ 
converges rather slowly to near the same stress of the tensile 
plateau. Again the material seems to undergo no Bauschinger 
effect. Fig. [6] re veals results of a similar experiment, conducted 
now with a bare copper film of 765 nm thick on polyimide, with 
no passivating coatings, strained at constant strain rate. Again 
o~l - o~t keeps constant for s > 0.01 and shows a more regular 
behaviour than the original data on cr L and crj. 

Fig. U\ displays and expanded view of the tensile branch ap- 
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Figure 4: Data of Gruber et al. (2) for a 251 nm thick copper film and d = 
297 nm on a polyimide substrate, strained at constant strain rate. The film has a 
passivating tantalum layer of 8 nm thick. The difference 07 - ctt between the 
two principal stresses keeps almost constant in the tensile region for e > 0.01. 
Upon compression, 07 - ctt go to a constant of the same magnitude. 
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Figure 6: Data of Gruber et al. (2) on a copper film of 765 nm thick and 
d = 423 nm on a polyimide substrate, strained at constant strain rate. The film 
has no passivating layers. The difference 07 - ctt of the principal stresses 
exhibits trends similar to those of Figs.[5]{5] 
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Figure 5: Data of Gruber et al. (2) on a copper film of 255 nm thick on a 
polyimide substrate, strained at constant strain rate. The film has passivating 
tantalum layers of 8 nm thick in both sides. The difference 07 - err between 
the two principal stresses exhibits the same trends of Figs.[3]and|4] 



pearing in Fig. [3] Together with the experimental data, the fig- 
ure shows the curve of true stress vs. engineering strain given 



by Eqs. ( 10 1 and (18} under the assumption that the engineering 
strain rate was kept constant in the tensile run. If the engineer- 
ing strain rate e eng is constant, the true strain rate s decreases 
with e eng , and the same should happen to the true stress cr. The 
curve was obtained from e eng = e eng /e - 1, where s - e VA is the 
function of cr x - cr y — cr L - crj given by Eqs. ( 10 1 and ((HI). Just 
to show an example of the procedure we adopted t c - 15 MPa 
and 2QT c /(7rds eng ) - 0.12. 

The two-dimensional plastic behaviour of submicron films 
has a special interest from a theoretical viewpoint. Making the 
same hypotheses and reasoning the same way as done to obtain 
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Figure 7: Expansion of the tensile branch of Fig.p] Circles represent 07, - o"t, 
with o"l and ctt taken from experiment. The solid line is given by Eqs. (To) 
and J8l with cr x - cr v = <tl - o"x and under the assumption that the engineering 
strain rate £ ena was kept constant all over the tensile run. 



Eq. ( [TO} , but now in a medium of three-dimensional polyhedra, 
one finds out that V • v + (5). Hence the three-dimensional 
model yields non conservation of the specific volume. As vol- 
ume variations can only be supported by the elasticity of the 
grains, one can recall Hooke's law and write V • v = -{1/B)p, 
where B is the bulk elastic modulus and /) the pressure vari- 
ation rate. This new equation determines that the stresses are 
not independent dynamical variables, plastic deformation is es- 
sentially a time-dependent problem and the stresses turn out 
dependent of the strain components [5J. Hence the three- 
dimensional version of the model yields strain strengthening 
or softening. The two-dimensional case is free of these com- 
plications, and thus sub-micron films are particularly suited for 
testing experimentally the hypotheses of the model for plastic 
deformation by grain sliding, which should hold the same in 



two- or three-dimensions. 
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